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ABSTRACT 


The interaction of two impinging plane turbulent jets has been 
investigated for the purpose of understanding the performance of propor- 
tional fluid amplifiers and momentum-exchange devices. Velocity and 
turbulence profiles were measured across and at various distances along 
the combined jet axis through the use of a constant temperature hot-wire 
anemometer. Experimentally determined jet deflection angles were found 
to agree quite well with those predicted in a previous work by the free- 
streamline theory. The results, which are presented in terms of normalized 
parameters, show that the jet interaction may be divided into three re- 
gions: mixing; transition; and fully developed. A vena-contracta ef- 
fect was found to exist in the mixing region. A discussion of each of 
the three regions is presented in detail along with generalizations for 


similar geometries. 
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I. INTRODUCTION 


In recent years, the need for a thorough understanding of the inter- 
action of two or three jets has become quite evident because of the in- 
creased application of jet interaction in fluid power technology. There 
are several important applications which may be cited herein. [In the 
field of fuel mixing, such as in jet or spacecraft engines, the intricate 
blending of two or more fuel jets is a crucial factor in determining the 
engine performance. Another example is found in the area of fluidics. 
The fluidic devices operate on a principle that one or more low energy 
control jets can impart sufficient deflection to a power jet and, thus, 
cause a significant change between the power levels of the receivers. 

The proportional fluid amplifier makes use of the most direct appli- 
cation of the interaction of three jets. Another fluidic element or 
logic device called "AND" gate operates through the interaction of two 
normally impinging jets. A simplified diagram of a proportional ampli- 
fier is shown in Fig. 1. The deflection of the main jet is achieved by 
an imbalance in the velocities or power levels of the two symmetrical 
control jets. Whenever a differential pressure exists between the control 
ports, the power jet is deflected by the momentum imparted by the control 
jets. This leads to a large or amplified pressure signal between the 
receiver or load ports. The interaction region of the amplifier is 
vented in order to avoid wall attachment due to the Coanda effect. Thus, 
the only external influence on the main jet is due to control flows. 

Although discovered in 1959, very little is actually known about the 
fluid dynamics of proportional amplifiers. A number of studies have been 
conducted on the overall performance of such amplifiers, but the general 


approach in these investigations hag been restricted to a "black box" 
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treatment rather than to a detailed understanding of the internal mechanism 
of the jet interaction. The realization for the need for a better under=- 
standing of the fluid dynamics of proportional or momentum-exchange de-~ 
vices has resulted, in the last few years, in the publication of several 
works (Refs. 1-4). 

Some basic research on jet interaction has been undertaken by Professor 
T. Sarpkaya in the Department of Mechanical Engineering of the Naval Post- 
graduate School. The present research is part of the investigation 
sponsored by the U. S. Army Research Office, Durham, North Carolina. 

The present thesis is a continuation of a project undertaken by N. 
Gungor (Ref. 5) to examine the interaction of two impinging jets. The 
evolution of the mixing region and the jet deflection angles were studied 
through the use of a constant temperature hot-wire anemometer. The veloc- 
ity and turbulence intensity profiles in the combined jet were obtained 


for various control and power jet velocity combinations. 
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II. EXPERIMENTAL EQUIPMENT 


The experimental apparatus together with the instrumentation used is 
shown in Fig. 2. The basic components of the equipment were: (a) jet 
assembly; (b) supply system; (c) hot-wire anemometer system; and (d) the 


hot-wire calibrator. 


A. JET ASSEMBLY 

The plexiglass jet assembly is shown in Fig. 3. In order to model 
part of a proportional amplifier, the axes of the two jets were made to 
intersect at 90 degrees. The jet-assembly design was essentially similar 
to that used by Gungor (5) with the exception of some additional modifi- 
cations to increase the accuracy and ease of measurements. 

The basic construction consisted of two 10 inch long, 1 inch deep, 
and 1/4 inch wide channels sandwiched between two 1/2 inch plexiglass 
plates. The control flow channel had a 1/4 inch setback as shown in Fig. 
6. The outer wall of each channel was cut to a sharp point (1/32 inch) 
to prevent attachment of the jet due to the Coanda effect. All sections 
were firmly assembled with dowels, and critical dimensions set with gage 
blocks to an accuracy of 0.001 inch. The channel surfaces were carefully 
hand-polished with rouge. 

The supply lines for each channel were 3/4 inch diameter Tygon tubing 
connected to special inlet devices to provide for the transition from cir- 
cular to rectangular flow channels. Wire screens were attached between 
the tubing and the inlet device by a flange connection. These screens 
served a dual purpose: (1) Eliminate any swirl that may have been im- 
parted to the flow in the supply tubing; and (2) Filter any foreign matter 


out of the flow. 
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The top and bottom plates had two slots cut out to enable the hot- 
wire anemometer mechanism to traverse along lines parallel to the combined 
jet axis. A telescope and slide mechanism were attached to the bottom 
plate. Degree marks were scribed (away from the interaction region) into 


the bottom plate, and a pointer was attached to the telescope slide. 


B. SUPPLY SYSTEM 

Air was supplied to the system from a 50 Hp. compressor at approxi- 
mately 200 psig. The air was first passed through a microfilter to remove 
oil, dust, and other impurities from the flow and then through the pri- 
mary pressure regulator (maximum input 400 psig and maximum output of 125 
psig.). Both the main and control jets were supplied air independently 
from this primary pressure regulator. 

Each jet supply line had a series of pressure regulators (10 to 130 
psig.) providing an accurately regulated flow into a Fisher-Porter rota- 
meter. Both the main and control rotameters were identical having a maxi- 
mum capacity of 14.6 scfm at 14.7 psia and 70°F with an FP=3/4-27-G-10 
tube and a 3/4-GNSVGT-54 float. Since actual flow rates were not measured 
at 14.7 psia, each rotameter was equipped with a pressure gage, and the 


flow rates were obtained by: 


Q = (rotameter reading (% of full flow)) x (14.6) x CF(P) 
where CF(P) = pressure correction factor obtained from a chart 


provided by the Fisher-Porter company for the rotameter. 


ce HOT-WIRE ANEMOMETER SYSTEM 
The Model 1050 constant temperature hot-wire anemometer manufactured 
by the Thermo-Systems, Inc. was used in the measurement of velocities and 


turbulence. The hot-wire probe used in the experiments was also manufactured 
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by the Thermo-Systems, Inc. (wire # 1272-T1.5). The wire had a cold 
resistance of 7.07 ohms and an operating resistance of 10.61 ohms. 

A special probe holder was constructed and is shown in Fig. 4. This 
holder fit into slots in the jet assembly which enabled the probe to move 
along lines parallel to the combined jet axis. A micrometer mechanism on 
the holder provided transverse movement of the hot-wire across the jet 


with an accuracy of 0.001 inch. 


D HOT-WIRE CALIBRATOR 

The hot-wire probe was calibrated with the Thermo-Systems, Inc. Model 
1125 Calibrator shown in Fig. 5. The calibrations were made with the 
high range calibration orifice having a diameter of 0.150 inch. The 
pressure in the calibration chamber was measured with water manometers 
with various ranges of pressure to increase the accuracy of the calibra- 
tion. The use of the calibrator will be discussed in the next section. 

A schematic drawing of the entire experimental system is shown in 


Fig. 7. 
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IIT, EXPERIMENTAL PROCEDURE 


The experimental procedure consisted of the following six main steps: 
(a) Calibration of the hot-wire probe; (b) Determination of the approximate 
jet deflection angles; (c) Adjusting of the hot-wire anemometer instrue- 
mentation; (d) Setting of the selected flow rates for each jet: (e) 
Measurement of local mean velocities and turbulence intensities; and (f) 


Data reduction. 


A. CALIBRATION 

The calibration consisted of the establishment of the exact relatione- 
ship between the anemometer voltage outputs and the corresponding veloc- 
ities. The Model 1125 Calibrator was used for this purpose. 

First, the hot-wire probe was placed about 1/16 inch above the ori- 
fice of the calibrator as recommended by the manufacturer (Ref. 6). Air 
was then let through the orifice, and the pressure in the orifice chamber 
was measured by a manometer. From these pressure readings, the velocity 


through the orifice was calculated by: 


[2P 
Un = Pp 


where P is the chamber pressure in psig., and fF is the air density. 
For each velocity the corresponding voltage output of the hot-wire anemo- 
meter was recorded. Finally, a calibration curve of velocity versus volts 
was plotted as shown in Fig. 8. 

The calibration curve in this form is not very useful since it cannot 
be used easily for a large number of measurements. In addition the de= 
rivatives of the calibration curve, which are needed for the turbulence 


intensity calculations, cannot be obtained from a calibration curve of 
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this type. A polynomial approximation is obviously necessary, especially 
for computerized data reduction. 

The hot-wire anemometer characteristics follow King's Law (Ref. 7) 
which is essentially a power law relation between voltage and velocity. 
A plot of log (velocity) versus log (voltage) data is shown in Fig. 9 and 
is seen to be very near to a straight line. A second degree polynomial 
approximation of this curve was found to fit almost exactly to the logari- 
thmic form of the data points. The derivation of this equation is given 
in Appendix A. The polynomial approximation to the calibration-curve is 
shown in Fig. 10 together with the actual data points for comparison. 

The calibration curve described above was checked frequently during 
the course of the experimentation. It was found that the wire performance 


remained essentially the same. 


B. APPROXIMATE JET DEFLECTION ANGLE 
In order to properly place the hot-wire probe near the center of the 
jet, the approximate jet deflection angle, oc , had to be determined. 
Two approximate methods were used by Gungor (5), but they were found to be 
rather time consuming and tedious. Since the accuracy required was only 
+ 2 degrees, a simple needle and thread method was used to determine oc . 
The apparatus consisted of a plastic holder with a long needle ate 
tached. At the end of the needle an eight inch piece of yarn was care= 
fully glued so as to have no preferred directional effects due to the 
method of attachment. The needle end was then placed at the center of 
the nozzle (Point M in Fig. 11). For each main and control flow combina- 
tion, the approximate jet angle could easily be read from the angle the 


thread made with the axis of the main jet. 


C. ADJUSTING OF ANEMOMETER INSTRUMENTS 

Prior to each experiment, the anemometer was turned on and allowed 
to warm up for about one-half hour. The cold resistance of the wire was 
recorded and checked to make sure that it had not deviated from previous 
values. Next the operating resistance was set into the device in accord=- 
ance with the operating instructions (Ref. 6). The system was then switched 
to the "RUN" operation mode. 

The fluctuating signal was monitored by an oscilloscope throughout 
the experiment. The internal feedback stability was adjusted by the "TRIM" 
dial on the anemometer to render the system stable within the range of 


velocities used. 


Dy SETTING OF FLOW RATES 

The channel flow rates were measured by the main and control rota- 
meters and adjusted by two pressure regulators. The main jet was main- 
tained at an average velocity of 87.5 ft/sec which corresponded to a 
rotameter setting of 60% of full flow and a Reynolds number of 11,140. 
The control jet flow rate was varied in the experiment to obtain differ- 
ent deflection angles. The control flow rates used and the corresponding 


velocities and Reynolds numbers are tabulated below. 
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Control Flow Rate Cc Reynolds 
__@ full flow) | a5 | ___ Number 





width of the main channel. 


E. MEASUREMENT OF VELOCITY AND TURBULENCE INTENSITIES 

Velocity and turbulence-intensity measurements were made along lines 
perpendicular to the jet axis at nine different distances from the inter- 
section of the inner walls of the jet channels. These nine stations were 
chosen to give both an accurate picture of the mixing region and the fully 
developed region. The following station distances were used: 0.15 inch; 
0.30 inch; 0.45 inch; 0.60 inch; 0.80 inch; 1.00 inch; 1.50 inchs; 2.50 
inches; and 4.00 inches. 

The probe was positioned using the telescope assembly discussed pre- 
viously. The telescope was first set to the approximate jet deflection 
angle. The probe was then centered at a particular station with the use 
of the telescope crosshairs. 

Once the probe was positioned, the traversing limits were determined 
to give an accurate picture of the jet profile. The hot-wire probe was 
then moved at suitable increments along the traversing line, and both the 
wire voltages and the R.M.S. voltages were recorded. 

A schematic drawing of the traversing line together with the appropri- 


ate nomenclature is shown in Fig. ll. 
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Fa DATA REDUCTION 

All data was reduced and plotted using the IBM-360 computer. The 
main data reduction program is presented in Computer Program B. Veloc= 
ities and turbulence intensities were calculated using the polynomial 
calibration-curve approximation previously described. 

The turbulence intensity was calculated by; 


= 


U dE 


where io is the R.M.S. voltage reading and dU is the derivative of 
dE 
the calibration curve evaluated at the value of the local mean velocity 


U. The derivation of the above equation is given in Appendix B. 
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IV. DISCUSSION OF RESULTS 


The main purpose of this investigation was to gain a better under- 
standing of the mixing of two turbulent jets. This interaction was found 
to be dominated by three regions: (1) the mixing region; (2) the transi- 
tion region; and (3) the fully developed region. The important character- 
istics of each of these regions will be discussed in terms of the velocity 
and turbulence intensity measurements. General similarity in jet inter- 
action for similar geometries will then be discussed through a comparison 
of the data presented herein with those obtained by Gungor (5). 

The most important characteristic in any jet interaction study is 
obviously the jet deflection angle. Figure 12 shows the comparison of 
experimental jet deflection angles with those predicted by the free- 
streamline theory (5). Also plotted on this graph are the jet deflection 


angles based on a simple momentum equation: 


U_\ 2 
ae 
| 
m 


It is apparent that the predictions of the free-streamline analysis 


<= pant 





and the experimental values agree quite satisfactorily. Deflection 
angles based on the momentum equation do not agree well at all. This is 
to be expected since the equation does not take into account the control 
channel setback. 

Gungor (5) also measured the deflection angles through the use of a 
model half the size of that used in the present investigation, i.e., 
W = 1/8 inch, aspect ration = 8. The table below gives a comparison of 
the two sets of deflection angles. Both sets of data have similar geome- 
tries and comparable Reynolds numbers with the exception of the aspect 
ratios. 
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Flow Ratio 
Q./Q. 





It is obvious that the deflection angles are almost identical al- 
though the actual magnitudes of the velocities were quite different in 
each case. Apparently, the deflection angle depends only on the velocity 
ratio rather than the actual magnitudes of the velocities. 

The three regions of the jet interaction can be seen clearly in Figs. 
13-21. In these graphs the velocity and turbulence intensity profiles 
for equal main and control jet velocities are plotted for various dis- 
tances along the jet axis. The mixing region is near the nozzle exit, 
and the two flows remain essentially separate with two velocity peaks. 

At larger distances from the exit, the two jets mix together until a 
single symmetrical jet emerges at a distance of about 6W as shown in Fig. 
19. This is the fully-developed region. The area between the mixing and 
the fully-developed regions is the transition region. 

Also shown in this set of graphs are the turbulence intensity pro- 
files. This combination of velocity and turbulence profiles shows graphi- 
cally the two turbulent mixing regions at the outer edges of the jet. 
These mixing regions are caused by the free shear layers. A third turbu- 
lent mixing region exists in the central core between the two jets and 


will be discussed later. 


Z2 


A complete set of data, for nine stations along the jet axis, is 
presented for U/U., = 1.00. Only two representative velocity and turbu- 
lence intensity profile combinations are presented for each of the remain- 
ing velocity ratios (Figs. 22-31). This is because all of the control 
flows exhibited similar characteristics. Turbulence intensities for 
U/U = 0.1605 have not been presented because of the inaccuracies in the 
derivatives of the calibration curve for low velocities. 

The mixing-region normalized-velocity profiles, for each flow ratio, 
are presented in Figs. 32-37. From these graphs the actual mixing pro- 
cess of the jets can be seen quite clearly. Close to the channel exit, 
the two velocity maximums are far apart with a large velocity minimum or 
"trough" between them. This preservation, for a short distance, of the 
identity of the main and control jets has been attributed by Douglas and 
Neeve (2) to a “brick wall effect". Further down the jet axis, the peaks 
move closer together as the momentum transfer across the turbulent shear 
layer decreases. Attemps were made to find some kind of similarity vari- 
able for the mixing-region velocity profiles but to no avail. Similarity 
will probably be obtained by analyzing each velocity mixing region sep- 
arately as mentioned in Pai (Ref. 8). 

In the mixing region the turbulence intensities are higher on the 
control jet side than on the main jet side. This higher turbulence in- 
tensity is caused by the low local mean velocities and not by the R.M.S. 
value of turbulence. With decreasing control flow, the difference be- 
tween the control and main jet turbulence intensity peaks rise which sup- 


ports the above reasoning. This same trend is also verified by Gungor 


(>) 
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As previously mentioned, there are three turbulent mixing regions. 
The graphs cited above presented only the turbulence intensities in the 
exterior regions of the velocity profiles. Figures 38-40 show, on an 
extended scale, the internal turbulence profiles for three velocity 
ratios. The intensity of internal turbulence decreases with increasing 
distance from the channels, as would be expected. By comparing Figs. 

38 and 40, it is seen that the greater the velocity difference between 
the two jets, the greater the internal turbulence intensity. This in- 
ternal turbulence is a result of the mixing of the two jets due to vor- 
tex motion. The turbulence decreases because of the vortex decay in 
the shear layer. 

It is apparent from the definition of turbulence intensity (see 
Appendix B) that it also represents the relative noise in the jet. 
Moynihan (Ref. 3) took the total pressure ratio as a qualitative measure 
of noise in a jet interaction. In graphs similar to those presented 
here, he showed three distinct noise levels in the jet of about equal 
magnitude. The results presented herein show that they are, in fact, 
not equal. The amount of turbulence ( aman is much less in the internal 
shear layer than in the external mixing regions. In proportional amplifi- 
er design, this fact is very important to keep in mind. 

The transition region has been defined as an area where the veloc- 
ity profiles do not exhibit dual velocity maximums nor are they symmetri- 
cal as in a fully-developed jet profile. A good example of this type of 
velocity profile is shown by Curve A in Fig. 43. Using crude terminology 
to explain the transition process, the jet "pops" from the mixing region, 
through the transition region, and into the fully developed region. The 


length of this transition region, for all control flows, was less than 
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one channel width. This region is short because of the characteristics 
of the vortex decay in the inner shear layer. As time increases, the 
turbulent vortices decay very quickly, causing a very short transition 
before a fully developed jet appears. 

In Figs. 41-46, similarity profiles are presented for the fully- 
developed velocity region. These graphs confirm the findings of Gungor 
(5) that fully developed profiles are similar when plotted against x/D. 
For convenience the transition region profiles have also been included 
in these plots. Since transition profiles are not symmetric, they can- 
not be expected to follow such a similarity trend. 

The effective pivot point of the jet, for each flow ratio, was 
examined by plotting jet maximum points along the jet axis similar to 
that shown in Fig. 11. Within the accuracy of this experiment, the pivot 
point was found to be approximately at the intersection of the centerlines 
of the main and control channels (Point M in Fig. 11) for all velocity 
ratios. For this reason, all of the deflection angles were measured 
relative to this point. This stationary pivot verifies the findings of 
Moynihan and Reilly (4). 

One of the most interesting plots of the mixing-region velocity data 
is shown in Fig. 47. This graph shows that the ratio of the two velocity 
maximums for a given flow ratio remains nearly a constant in the mixing 
region. This fact could be very important in attempts to modify the free- 
streamline theory to better describe the mixing region. It is also seen 
from these curves, that the length of the mixing region increases with 
increasing control flow. The maximum length is 6W for Uele= 1.00. 

The dashed lines to the right of each curve represent the transition 


region and show how small this region actually is (about 0.5W). 
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Figure 48 shows the acceleration of flow and the presence of a vena- 
contracta in the mixing region. Here the ratio of the largest maximum 
velocity in a particular profile to the maximum velocity in the main jet 
channel has been plotted against D/W. A maximum occurs at a distance of 
2.5W for almost all flow ratios. This may be interpreted as a "pinching" 
or vena-contracta of the combined jet. 

This vena-contracta effect is caused by two phenomena. First, the 
control jet "pinches" the main jet because of the direction of the initial 
control jet momentum. Secondly, as the jet moves further from the nozzle, 
surrounding air is entrained by the turbulent mixing. Eventually, this 
entrainment and mixing causes the jet to expand and decrease in velocity. 
The net result is a vena=-contracta type of velocity profile. The amount 
of pinching of the main jet is certainly related to the relative momentum 
of the control flow. The vena-contracta is almost non-existant for low 
control flows. 

One of the primary purposes of this investigation was to determine 
if any general similarity in jet interaction existed for similar geometries. 
In order to accomplish this, the results for W = 1/4 inch were compared, 
wherever possible, with the results of Gungor (5) for W = 1/8 inch. 

A comparison of centerline velocities, for various distances along 
the combined jet, is presented in Figs. 49-54. The paca pl ratio is 
seen to have the same general functional relationship with D/W even though 
velocities and channel widths are quite different. Figure 55 is a graph 
of D/W (corresponding to the maximum centerline velocity) versus the flow 
ratio. It appears that the centerline velocity is a function of D/W and 


U/U, for similar geometries and comparable Reynolds numbers. 
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V. CONCLUSIONS 


The following conclusions can be drawn from the data presented here- 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Experimental jet deflection angles agree with the predictions 
of the free-streamline theory. 

The jet deflection angle is independent of the magnitudes 
of the main and control jet velocities for similar geo- 
metries and depends only on velocity ratio. 

For two impinging jets, there are three turbulent mixing 
regions with the turbulence intensity directly dependent 
on velocity differences. 

Fully-developed jet profiles for each flow ratio are 
similar when plotted against x/D. 

The effective pivot point of the deflected jet is the same 
for all flow ratios, i.e., at the intersection of the 
channel axes. 

In the mixing region, the ratio of the two velocity maxi- 
mums is nearly constant for each flow ratio. 

The length of the mixing region increases with increasing 
control flow. 

A vena-contracta effect is present in the combined jet at 
distances of about 2.5W from the channel exits. 

The centerline velocity in the combined jet is a function 
of only D/W and velocity ratio for similar geometries and 


comparable Reynolds numbers. 


2/ 


VI. 


The investigation described herein should be extended to the study 


of the following areas: 


(1) 


(2) 


(3) 


(4) 


(5) 


The effect of different setbacks and aspect ratios on 

the jet deflection, velocity, and turbulence distributions. 
Effect of low Reynolds number or laminar flow on the jet 
interaction. 

A more thorough investigation of the turbulence in the 
jets using two-wire correlation to check for dominant 
noise frequencies and self-preserving turbulence. 
Verification of theoretical predictions for similarity in 
the three velocity mixing regions. 

An extension of the free-streamline theory approximation 
for the mixing region of two impinging jets to include 


two unequal free stream velocities of constant ratio. 
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APPENDIX A 


LOGARITHMIC CALIBRATION CURVE APPROXIMATION 
The basic relationship between the voltage and velocity for a hot- 


wire anemometer is given by King's Law: 


(1) <AE=«kt" 


where 
AE = differential voltage 
U = local mean velocity 
K,n constants 
taking the logarithm of (1) 
(2) in( AE) = 1n(K) + nin(0) 


letting 


y = In( AE) 


iH 


V = 1n(0) 


and approximating (2) with a polynomial one has 


2 m 
(3) V=a, + ay + ayy = onesies Aer ae 


Various values of m were tried using a Chebyschev Curve-Fit as shown 
in Computer Program A. A second degree (m = 2) polynomial approximation 
was found to give the best fit to the actual data (see Fig. 10). 


The coefficients a for the polynomial were found to be: 


a 

m 
a, = 3.494994 
ie 3.008039 
a, = -0,003591135 
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APPENDIX B 


DERIVATION OF TURBULENCE-INTENSITY EQUATION 
For a turbulent flow, the instantaneous value of the anemometer out- 


put is given by: 


(1) E=E+Ve~ 
where 
E = voltage output corresponding to the mean velocity 
at a given point (U). 
a = R.M.S. value of the voltage output corresponding 


to the R.M.S. value of the fluctuating component 
of velocity ( ywuws). 


Thus from equation (1) we can write the equation for the instantaneous 


velocity 


(2) vu=T WT? 


and 
(3) E=f(U) = f@U + ra 
also 


(4) dE = £'(U)dU 
since the fluctuating component of velocity, re is relatively small, 
we can expand equation (3) in a Taylor Series neglecting higher order 
terms, 

(5S) £f¢ g Az) = £(0) ra Preesareiy. + eee 

therefore 

(6) E eE+Ve" = £(U) alte 

by definition E = f(U) so that 


7) Ye? lee ® 


and 
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(3) Ware es Ven 


£' (0) 
the turbulence intensity is defined as 


(oye tem mie. 


U U £°(0) 


since £'(U) = dE 


d 


the final expression for turbulence becomes 


U U dE 








Equation (10) is easily adaptable for computerized data reduction 


as shown in Computer Program B. The values of U, and dU/dE are obtained 


from the polynomial calibration curve. 
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